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i—i \ Abstract 

In the minimal SU(5) SUSY GUT d = 5 operators lead to p — > K + u decay 
X/} \ with the proton life time of the order of 10 28 years for the natural choice of the 

■ parameters of the theory. This value is in strong contradiction with experimental 

t— I . bound Tp^,Ku > 10 32 years, d = 5 operators are induced by colored Higgsino 

exchanges and are closely (through SU(5) and super symmetry) related to another 
wrong prediction of SU(5) SUSY GUT: md/m s = me/m^. We demonstrate how in 
the model where reasonable pattern of quark and lepton masses and CKM mixing 
angles are obtained proton decay can be suppressed and proton life time can be 
' close to the present experimental bound. 

o 

00 ■ 

1 Introduction 

P-T 

As everyone is aware of the most attractive candidate for the physics beyond the Standard 
Model is low-energy supersymmetry (see reviews 0). It helps in solving hierarchy problem 
- so, GUT's get firm theoretical foundation. Experimental signal in favor of SUSY GUT 
>• ■ comes from the numerical value of electroweak mixing angle, sin 2 6 exp ~ 0.23, which 
nicely coincides with SUSY GUT value, while contradicts non-SUSY GUT prediction 
sin 2 6 « 0.21 ||. Another manifestation of this phenomena is prediction of the a s [Mz) 
value which nicely coincides with LEP and other low-energy measurements. 

One of the most spectacular prediction of Grand Unification is proton decay. In 
nonsupersymmetric theories proton decay through d = 6 operators mostly via p — > e + 7i° 
channel. Modern experimental bound on this particular mode is Tp_> e + w o > 10 32 years 
which strongly contradicts prediction of the SU(5) GUT: r p _^ e + n o = 10 28±2 years [|]. In 
super symmetric GUT's operators with d = 6 are also generated. But since in SUSY GUT's 
unification scale is approximately 30 times larger than in non-SUSY GUT's, proton life 
time due to operators with d = 6 is of the order of 10 34 years []J which is beyond discovery 
possibilities because of background problems. 

In supersymmetric models the d = 4 trilinear B and L violating couplings can be 
introduced. They mediate fast decay of proton. That is why one should impose on the 
theory condition of absence of such operators. 

However, the d = 5 operators are induced in SUSY GUTs by exchange of the color- 
triplet Higgsinos, which are partners of the Higgs doublets in the GUT multiplets [|J. 

x Talk presented by M.Vysotsky at the Conference "Quarks 98", Suzdal, Russia, May 18-24, 1998. 
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In the second section of this paper old result for r p _> K + u in SU(5) SUSY GUT 
will be reanalyzed. Feynman diagrams which induce this decay are shown on Fig. 1. 
This refreshment is necessary since at the time when || was written on the one hand 
lower experimental bound on t p ^k+ u was two orders of magnitude weaker than now, 
on the other hand, neither m top was known nor Kobayashi-Maskawa mixing angles were 
measured with modern accuracy. Using these updated numbers together with old (but 
still valid) value of decay matrix element from || we get our central statement: proton 
decays too fast and we cannot naturally be within experimental bound 0: r p _ iK + L/ > 10 32 
years. In this way we come to the following conclusion: minimal supersymmetric SU(5) 
GUT should be modified. In other words, some mechanism for suppression of the nucleon 
decay amplitude is necessary and we are not the only who shared this point of view . 

In SO (10) theory suppression of the proton decay can occur for the following two 
reasons: 

1) The scalar sector is arranged in such a way, that the nucleon decay parameter - 
(M^ 1 )n vanish or is strongly suppressed |§. SO(10) model in which the proton decay 
is strongly suppressed and in which gauge and fermion mass hierarchy was explained 
naturally was suggested in |J. 

2) Another possibility of stabilizing the proton by implementing the 45-plet with VEV 
towards Tr direction in the Yukawa sector was suggested in refs |L0| f\ and fTT|j . 

Our approach here is different. We study SU(5) SUSY GUT where unsatisfactory 
relations of minimal theory = m s and m e = at GUT scale are avoided. 

These unsatisfactory predictions come from the same Higgs-matter multiplets cou- 
plings which generate d = 5 operators. Having in mind valuable way to solve the mass 
degeneracy problem we will work on d = 5 operator generated proton decay in this scheme. 
One can see that what is done in the present paper for the suppression of the nucleon 
decay in the framework of SU(5) theory is analogous to what was proposed for SO(10) 
SUSY GUT in |TD[ and [JTTJ . It is natural to assume, that the renormalizable couplings 



with Higgs multiplets has only the third generation, while the lighter generations get 
masses through higher order terms, by the mixing with there nearest heavy neighbours 



13]; we suppose that this higher order terms for up quarks are antisymmetric. Since qqT 
coupling is symmetric in the generation indices it vanishes for the light generations and 
exists only for the third family. This leads to the strong suppression of the nucleon decay 
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The paper is organized as follows: In part 2 we present the d = 5 operators for general 
SU(5) theory. Part 3 deals with proton decay in minimal SUSY SU(5). Part 4 contains 
the solution of m ; = m d problem in our extended S77(5) SUSY GUTs. In part 5 we 
consider proton decay in our model and part 6 contains discussions and conclusions. 

2 d = 5 Operators 

Fermion sector of the SU(5) SUSY GUT consists of the one pair of fermion supermultiplets 
5 + 10 per generation: 

5 a = (d c , l) a , 10 = (u c , q, e% , (1) 
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The analyses of the nucleon decay in this model was presented in [fL2| 
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where a = 1, 2, 3 is a family index. 

The Higgs sector contains the following chiral supermultiplets: E ~ 24 in adjoint 
representation of SU(5) and 5 and 5-plets H, H: 

F = (T, H u ), H=(T, H d ) . (2) 

The SU(5) invariant Yukawa couplings which generate masses of the up and down 
quarks and charged leptons are respectively: 

10 • f u ■ 10H, 10 • f d ■ 5H, (3) 

where T u and T d are Yukawa coupling constants (family and SU(5) indices are suppressed). 
Decomposition of these couplings in general have the form: 

r u 10 ■ 10 ■ H qY u u c H u + qAqT + u c Be c T , (4) 

r d 10 • 5 • H -> qY d d c H d + qClf + u c Dd c T + e c YJH d . (5) 

After integrating out the colour Higgses T, T with masses of the order of Mqut we 
obtain the following d — 5 operators : 

L = -^—(qAq)(qCl) , (6) 

R = J ^—(u c Be c )(u c Dd c ) . (7) 

In general r U)< j (see (j^)) can be some functions of E/M, where E breaks SU(5) down 
to the G 123 = u\l) Y x SU(2) W x 5£/(3) c group at the scale M GUT ~ 10 16 GeV. M is 
some fundamental scale, M 3> M GUT . 

After diagonalization of the Yukawa matrices by biunitary transformations: 

L+Y U R U = Y^ iag , L + d Y d R d = Y d Dtas \ R + e Y e L e = Y e Dia ° (8) 

and by proper redefinition of the quark and lepton fields all operators can be rewritten 
in the mass eigenstate basis, where the interaction of the quark-lepton superfields with 
colour-triplets have the form: 

qAqT - q u L+AL* d q d T , (9) 

u c Be c T -> u c RlBR e e c T , (10) 

qCIT -> q u LlCL e l e T - q d L+CL e l„T , (11) 

M c /)d c f -> u c R T u DR d d c T ; (12) 
in this basis the current-gauge superfield interactions have the form: 

92qtVq d W^ +g 2 lthW^ (13) 

where V is the CKM matrix: 

V = L T u L d (14) 
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The Ol type d = 5 operators which induce proton decay with neutrino emission by 
exchange of wino lead to the following baryon number violating four-fermion interactions 
(we omit charge conjugation matrix in fermion braces throughout this paper): 

\r C S dv) ■ {<dt)i^ P )e ahc (15) 



M H3 



where 



C = C(I) + C(II) + C(III) + C(IV) , (16) 

C^\l) = -g!(LjCL e ), p (LtAL^V, a (V + Ul(u^ d") , (17) 

C^ du) (H) = -gl{L+AL* u ) lP V Pa {L T e C T L* d ) prT (V + Ul(u?, d°) , (18) 

C^iHI) = gl(L+AL d ) 5a (L+CL e ) 0p V Pl I(u^ &) , (19) 

Cjr a T\lV) = -9l{LtALlUV Pl {LtCL e ) 5p I{u^ e>) , (20) 
while the operators with charged leptons are: 

C ap8p Ue) ■ {K d b)( u l eP ) £abc (21) 



Mh 3 



where 



C = C(I) + C{II) + C(III) + C{IV) , (22) 

cS$™(i) = -gl{LZAL* d U{LiCL e ) ap <y+ui(fr,*) , (23) 

C%% ue) (II) = 9 mAL d U(V + ML^CL e ), p I(d^ up) , (24) 
C$™(M) = ^(L+(7L e ) 5p (L+iL*) 7CT ^(^ + ) CTa /(^,^) , (25) 

c%$!r\iv) = gmAL:) Sl (v + ) ja (& T L:) p ^i(d\ (26) 

(a, and a,b,... are the family and colour indices respectively). C(J) and C(III) 

correspond to the vertex diagrams (see Figs. 2a,c,e,g) while C(II) and C(IV) to the box 
diagrams(see Figs. 2b,d,f,h). In (|l^)-(^) I denotes the result of the integral over the 
loop and is given by the following formula [pOfl : 



'(u, i) = T^^^ [-^r I- 4 - 4) . (27) 

167r z m~ — m T V m~ — mfj. md, m ~ — mf?. mf ? . / 
and analogously for I(d, u) and I(u, e). 
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For degenerate squark masses from (|27|) we get: 



where 



1 Jx~Z 1 

I(u,d) -> I(m^,x w ) = — — — ( Xw In x w - x w + 1) , (28) 

lo7r^ (1 — x m y 



*- = (^j • (29) 



The function I(rriq,x w ) has the following behaviour: 



In i„ 



if 2T ro > 1 



167T 2 mqI(mq,x w ) = \ 0.5 if = 1 (30) 

\J Xw if < ^ 1 

In order to estimate stop contribution the following relations will be used: 

m 2 = m~ + m 2 , (31) 



T(i 3\ = 1 m w m\x w \iax w + (m| + m 2 )(l - a; w ) ln(l + m 2 /m 2 ~) 
l ' j 16vr 2 m 2 m 2 (l-x 5i ,) + (l-x w ) 2 m? 



3 Proton decay in minimal SUSY SU(5) 

In minimal SU(5) T u and are SU(5) singlets and the following relation between the 
Yukawa matrices occurs at the GUT scale: 

Y U = A = B , Y d = Y e = C = D . (33) 

In addition in the minimal SU(5) theory f ^ = f u which leads to equality of matrices 
which are used to transform u and u c to mass eigenstate basis. Therefore in the mass 
eigenstate basis these couplings have the form: 

r u 10 • 10 ■ H -> qY® ia9 u c H u + q u Y u Dia9 Vq d T + u c Y° ia9 V*e c T , (34) 

T d 10 ■ 5 ■ H -> qY d Dia9 d c H d + e c Y e Dia9 lH d + q u V*Y d ia9 l e f- 

q d Y^ ia9 l v f + u c V*Y d Dia9 d c f , (35) 

where V is the CKM matrix. 

Taking into account the renormalization effects between the GUT scale and the SUSY 
scale the equalities ( |3"B| ) are violated. Using the ([53|) as a boundary conditions, at the 
SUSY breaking scale the A <g> C and the B ® D from flf]) can be expressed by the Y u ® Y d 
product: 

A®C= {A S ) L {% ® Y d ) , (36) 
B®D = {A s ) R {Y u ®Y d ) . (37) 
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The As coefficients describe the renormalization effect between the GUT and SUSY break- 
ing scales. The numerical effect of As factor will be discussed in what follows. 

Calculation of the amplitude of proton decay consists of two steps: calculation of 
1-loop Feynman diagram(s) where transition between sparticles and particles occur and 
calculation of matrix element of the corresponding four-fermion operator between proton 
and Kv system. Quarks and leptons of the second and/or the third generations give the 
main contribution. As these particles (except neutrinos and s-quark) do not participate 
in proton decay, their scalar superpartners go into the loop and are transformed into light 
species through wino exchange. Diagrams which describe the decay p — ■> Kv^ are shown 
on Fig. 1. Sum of the amplitudes which are given by diagrams shown on Fig. lb and 
Fig. lc equals zero, and we are left with four diagrams shown on Fig. la. 

As it usually occurs in SUSY models, vertices on the diagrams shown on Fig. 1 
are known - they are the same as in nonsupersymmetric SU(5). Less is known about 
propagators - masses of squarks and wino. Let us remind, that two Weyl higgsinos from 
two Higgs doublets with unit charge mixes with two Weyl winos and two massive Dirac 
particles are formed. Mixing matrix contain four parameters /i, Myy, 5 , 2^i/2, 92^2/^ which 
are constrained by one equation: 

{ my + { m ? = Mir . (38) 

In SUSY GUT = a 2 /a 3 m g , where trig is gluino mass. Since from Tevatron bounds 
gluino should weight several hundreds GeV at least, we have m-^ ^> g 2 Vi/2, #2^2/2 and 
instead of dealing with two mass eigenstates in box diagram we could take into account 
only jy-bosino exchange. 

Calculating diagrams of Fig. la we obtain 

M = 2 L usyy2J Uu d)AaA L x 

( Vl /V2) ■ (V2/V2) Mjj 3 

[{v L d a L ) (u b L s c L ) + (u L s a L ) (u b L d c L )]e abc , (39) 

where for triplet higgsino-matter coupling constant we use / = m q Vik/(v/y/2). W-bosino 
transform sparticles into particles with the constant g 2 - Factors A take into account short 
(As) and long (Al) distance renormalizations of decay amplitude. The index I of A l s 
refer to the contribution of the two light generation particles propagating inside the box 
diagram of Fig. la . Factor Al is the long-range renormalization factor due to the QCD 



interaction between the SUSY breaking scale and 1 GeV scale [19 



Al = ( a 3 (lGeV) \ ~ 2/3 /«3K) \ " 18/25 f a 3 (m») \ ~ 18/23 
L \ a 3 (m c ) J \a 3 (m b ) J \a 3 (m z ) J " 

and for a 3 (mz) = 0.120 using a 3 running at two loops we get Al = 0.32. 

In (|39D factor 2 comes from 2 diagrams, are the elements of Kobayashi-Maskawa 
matrix, v\ and v 2 are Higgs fields vacuum expectation values and is mass of Higgs 
triplet, M kz w M GUT = 10 16 GeV (let us remind that for minimal SUSY SU(5) Higgs 
triplets interactions with quarks and leptons which generate operator Ol are described 
by the Kobayashi-Maskawa matrix). 
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It is convenient to rewrite (|39"D introducing angle /3, tan/5 = vi/v 2 and expressing 
■u 2 + v\ through G F : 

M = 8 ^ G M™2m dVUS ^ A ^ L dl) (u b L sl) + (v L sl)(u b L dl)}e abc . (41) 

For the matrix element of operator fl4l|) between hadronic states we use the result 
obtained in : 

(vK+\e ab °[(v L dl)(u b L sl) + {v L sl){u h L dl)]\p) = f ^^,,J uP L p)K , (42) 



(M A + M s )/2 



where P L = ~(1 + 75). 

Finally, from (0) and © we get: 



^ 16 • f3GG F m c m s g 2 2 V cd V us , 

p ~* Ku = M^sin2/5[(M A + M E )/2] /(m ^ ^)^(^) = *(^*p)* ■ (43) 

Short distance renormalization factor As depends on the numerical value of tan/5 
20f . Making an attempt to suppress the proton decay amplitude we take the value of 



A 1 

tan/5 which minimize the ratio sin |g ; so we use sin/5 = 0.965 (tan/5 = 3.68), A s = 1.4, 

A 1 

|-3 = 2.79 (detailed calculations will be published in an extended paper). 



sin 2/3 

Substituting numbers in (EST) we get 



(ml - m\f 2 _ / m c m s 10 16 GeV /5 

32vrm3 X ~~ { 1.3 GeV 175 MeV M kz 0.007GeV 3 X 

I(mg,x w ) A L A l s 0.51 ^ 1 



J(500GeV, 1) 0.32 1.4 sin 2(3 ' 4.5 • 10 27 years 



(44) 



Modern experimental bound is t p ^ku > 10 32 years 0. Variation of parameters could 
hardly help in enhancing proton lifetime that much, q and W with mass scale several 
TeV did not seem appealing, neither is w 10 18 GeV (let us remind that SUSY GUT 
unification scale is M GUT = 10 16 GeV and M^ 3 w X/ g ■ M GUT , where g is gauge coupling 
at unification scale, while A is a constant of Higgs multiplets selfinteraction) . 

The proton lifetime for the different values of squark and wino masses are presented 
in Table |l] . As we see the reasonable lifetime is obtained for m§ = 5 — 10 TeV and m^ 
about 1 TeV or less ( in this domain x w is small and the function mgl(mq,x w ) can be 
described by the asymptotic formula ([30|)). If one wants to have lighter quarkino, with 
mass less then, say 1 TeV, then proton decay should be somehow suppressed. 

Before we will go to the main part of our paper let us estimate how much the contri- 
bution of the third generation particles in the proton decay amplitude can be. If instead 
of §l (P>l) on Fig- la we substitute bi (tl), we will get the following extra factor in the 
amplitude (|39|): 

A h s (b)m b .V ub . 4.1 -4.5 0.002-0.005 

4 , — — = x = 0.1 — 1 . (45) 

A l s m s W U s 0.1 -0.3 0.22 v ; 
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Table 1: The value of proton life time in Standard SUSY SU(5) in units of 4.5 -10 27 yeares. 
Allowed domain of m^, rriq values is in low left corner. 



fflq 


100 GeV 


200 GeV 


500 GeV 


1 TeV 


5 TeV 


10 TeV 


100 GeV 


0.04 


0.03 


0.04 


0.07 


0.54 


1.48 


200 GeV 


0.31 


0.16 


0.13 


0.17 


0.84 


2.14 


500 GeV 


7.7 


2.6 


1 


0.78 


1.8 


4.0 


1 TeV 


1.1 • 10 2 


30.7 


7.8 


4 


4.2 


7.4 


5 TeV 


6.3- 10 4 


1.6 • 10 4 


2.7- 10 3 


7.7- 10 2 


10 2 


78 


10 TeV 


10 6 


2.5 • 10 5 


4.1 • 10 4 


1.1 • 10 4 


7.8 ■ 10 2 


4- 10 2 
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Ag is the short range renormalization factor for the heavy generations and Ag(b) = A l s . 
Stop substituted instead of cl on the upper line of Fig. la lead to the following factor 
in the amplitude (|39|): 

A h s (t) I(t,d) Vtm t .V td V tSl 



A l s I(u,d) m c V cd 
1.9 6.1 • 1(T 6 2.4 ■ 180 (0.004 0.014) x (0.034 0.046) 



0.3 -=-1.5, (46) 



1.4 6.3 • 10- 6 1 -=- 1.6 0.22 

where we use nig = = 500 GeV. Ag(t) = 1.9. 

Let us stress that amplitude ( ffl]) is defined at \i = 1 GeV. Since t-quark mass is not 
renormalized from the virtuality which equals to its pole value mt = 180 GeV to virtuality 
1 GeV a compensation factor r\ t should be introduced in A L : 



m t (l GeV) 

Vt = 

m t 



a 3 {l GeV) 



a 3 (m t ) 



(47) 



for a 3 (M z ) = 0.120 this factor equals 2.4. 

From (|45"|) and (|46|) we see, that for the maximum mixing between first and third 
generation allowed experimentally contribution of third generation particles into proton 
decay can compete with that of second generation. Amplitude with intermediate stop 
interfere with that with intermediate scalar charm quark and may partly cancel it; however 
compensation with 1% accuracy which is needed to satisfy experimental bound t p ^„k > 
10 32 years is unnatural. 



4 Predictive ansatz for Yukawa couplings and sup- 
pression of proton decay 

By focusing on the fermion mass pattern, it is natural to suggest that only the third, 
heaviest family acquires masses through ordinary renormalizable Yukawa couplings, while 
the mass terms of other families appear from higher order (may be Planck scale) operators, 
which can be generated by heavy particle exchange mechanism [tH] : 

= -C ■ 10 ■ im + -— ■ 10 • 10 • EH + ^- ■ 10 • 10 ■ T?H , (48) 
Y 4 AM M 2 K ' 

b' - a' 

W$ = V2d ■ 10 ■ 5H + V2— ■ 10 • 5 • £ • H + • 10 ■ 5 • S 2 ■ H , (49) 

where C , B',... are matrices in generation space. 

In order to be closer to the realistic mass matrices let us suggest for them the following 
form: 

C'aP ~ ^3a^3/3 , B' ap ~ 5 3a 5 2f S + k B 6 2a S 3 /3 , A' af3 ~ + k A 5 la 5 2 p , 

d a p ~ 5 3a 5 3/3 , b' a p ~ 5 2Q 5 2 f3 , a' a/3 ~ 5 2a 5i/3 + k a 5 la 5 2l3 . (50) 
In ( p0|) k are Clebsch factors. For Wy we get: 

10 x 10 = 5 + 45 + 50 , (51) 



9 



while for W Y we have: 

10 x 5 = 5 + 45 . 
For the bilinear Higgs fields product we have: 

24 x 5 = 5 + 45 + 70 



(52) 



(53) 



so the Ei7 could belong to 5 or 45 and in these cases B' is symmetric (ks = 1) or anti- 
symmetric {k,B = -1), respectively. However, because in 24 x 24 x 5 several invariants of 5 
and 45 plets and also 50-plet do occur, there exist many invariants and many possibilities 
for kA and fc a 's. 

In what follows matrices A', B' will be taken antisymmetric (kA = ks = — 1) and this 
will be crucial for the proton stability. 

In other words we suppose that for some reason the composite operators Eif and T?H 
are participate in expression (]4Sf) only in representation 45 . 

Insertion of E in higher order terms helps to avoid the degeneracy of the masses of down 
quarks and charged leptons. As it was assumed matrices B' and A' are antisymmetric, 
while a' is symmetric with respect to the family indices. (This can be attributed to some 
symmetry reasons). Then after substituting the VEVs of E, H and H Yukawa matrices 
for up and down quarks and leptons will have the forms: 



(54) 



Because according to our choice B' and A 1 matrices are antisymmetric in the family 
space while the qqT coupling is symmetric on the generation indices only 33 element of 
the matrix A is nonvanishing (as we will see this fact is crucial for proton decay): 





( 





.4 






f 





a>i 






( o 


0-2 




Yu = 




-A 





?)■ 


Y d = 




a i 




l \ 


*- 


a 2 


b 2 


!) 




\ 





-B 






{ 








d J 




K o 








A 



C 



( 





\ 















\o 





C ) 







a i 






a, i 


bj 


"") 










d y 






(55) 



(56) 



[|56|) depend on the SU(5) repre- 
belong. In numerical estimates 



The values of matrix elements of matrices (|5^ 
sentations to which higher order Higgs terms in I 
we will take di = a\\ concerning b\ two possibilities will be considered (see later). 

Structure of the matrices ( [5*4] ) resembles the ansatz proposed by Georgi and Jarlskog 
Lately a number of authors 



in an SU(5) GUT |22 



M p5[ have reexamined this texture 



in a supersymmetric context. From (E3) it is easy to find, that 



A w V ^uK , B w J \ c \ t , C w A t , a x w J \ d X s , 6i w A s , 



a 2 w J AeA^ , 6 2 ~ A^ , d = X b = A T . 



(57) 
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The Yukawa matrices are diagonalized by the transformations given in 
from fl54|) we have: 



where 

r(l2) 



/ cos9 L 



12 



sin 0fy 

V o 



Ll 23 ) • Li 12 ) , L n 



12 



— sin 9 

cos 0f 2 " 



'(12) 



r(23) 



L P = Li 12 ) 



/ 1 


V sin# : 



cos # 



23 



23 





— sin 6 

cos 9 



23 



23 



where for Y 



(5f 



(59) 



■(12) 



cos 12 
sin #f 2 d 




-sin#f 2 d 
cos^ 2 d 
1 



sin #f 2 " ~ 



AC 



sin^ ~ 



Therefore at the GUT scale we have: 



\ 



1 £ 

1 C 



1 

IP" 

B C 



L r , 



while for -R n! d je matrices we have: 

i a o 

p _ I __AC i _B 

-^tt — | £J2 X 



£1 £. 1 



L (12) 



( cos 6»f 2 e 
sin 6>f 2 e 
V 



cos# L 



15 

'e 
12 



sin #f 2 e 




1 



5 

C ' 



sin# 12 d ~ 



ai 

&7 ' 



^, sin^-? 2 - 



1 


«i 
6i 


^ 




( 


1 


012 

b-z 





ai 


1 





, L e — 




02 


1 











1 J 




\ 








1 



In 





1 
1 



/- 1 f 

02 1 

V 





1 



After diagonalization the Yukawa matrices have the form: 

Y U D ^ = (A u , Ac, A t ) , Y d Dia9 = (X d , A s , A 6 ) , y® ia Q = (A 



1 ) -^U) -^r y 



From (111), (0) and (1621) one can find the CKM matrix elements: 



(60) 



(61) 



(62) 



(63) 



(64) 



V u . 



cb 



V, 



lib 



v t 



ts 



v td = -< 



AtA s 



(65) 



rCXp 

'cb 



As we see on the GUT scale the value of the V co element is too large (V^ 
0.036 -i- 0.046). It appears [2"3fl-||25||, that the desirable relations between masses and 
mixing angles are obtained on the electroweak breaking scale after taking into account 
the renormalization effects. 
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5 Proton decay in extended SUSY SU(5) 



Let us estimate now the proton decay probability in our model. Let us start with p — > 
K + u, mode, which dominates in the minimal SU(5). This decay is described by the 
diagrams Fig. 2a-d. 

From ([H2]), ( |5H| ) it is easy to see, that (L^ AL* u ) a p exactly vanish for a = 1, 2. Therefore, 
the amplitudes (fHf), ([19|) and ( j20"|) do not lead to the proton decay as they produce 6-quark 
in final state. 

The amplitude described by eq. (|T7|) is suppressed for another reason: as we see from 
(|j~7l ) this amplitude do not vanish if a — 3. However in the inner line of diagram Fig. 2a u, 
c and t squarks run. Assuming for a moment that integral / in (|TTD is family independent, 



taking sum over (3 and using (|14|) we see that (|17D is proportional to (L^AL^) a3 and in 
the external line still the b quark is produceed. In this way we see that d or s quarks 
which can participate in proton decay are not emitted. However the above argument is 
valid only if the equality m„ = mg = m 4 ~ holds; heaviness of the top quark breaks last 
equality, so p — > Ku decay through diagram Fig. 2a do occur. So, taking into account 
the shift of / function the nucleon decay will take place due to heavy stop exchange, but 
the suppression factor ^j- will appear, where AI is: 



These arguments work for both p 
Introducing the parameter 



AI = I(t,b) - I(u p ,b), 
Kv and p 



(3 = 1,2 
tiv decays. 



(66) 



for the p 



Sg 2 2 GG F A h s (t)Ajl 
M^ 3 sin2/? 

Kv a decay widths f\ we get: 

(m 2 p - m 2 K f 



(67) 



r(p - Kv, 



i 

3271777.3 



-Ft 



m t \V ub \ x 



Vts(C*12 + ^22-^^12 

r(p 



2/3a 



(M A + M s )/2 



1 - 2/3a 



mi 



m 2 K f 



327rrrtp 



M A 



m 2 t \V ub \ 2 x 



(68) 



Vts(Cl2(L e ) 2 l + C*2l(Iv^) 1 2 + C22(-^!)l2(-^e)2l) 



2/3a 



(M A + M s )/2 



while the widths of the p 
T(p 



Vtd(C 2 l + C 2 2(-^e)21 

• nu a decays are: 



1 - 2/3a 



M A 



m r . 



it v„ 



32nm 2 1 



A/\ 2 
/ ) 



V2_ 



m 2 \V ub V t d\ 



Cyi + C22{L~ d ~ )l2 



(69) 



(70) 



'necessary matrix elements are presented in the end of this section 
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Table 2: An order of magnitude estimates of the proton partial life times in the units of 
10 32 years. 



C 22 = 



C 22 — 2A. S 



V - 

V - 
V 
V 



Kv e 

• 7TZ/ fi 



' 10 
200 

> 10 
300 



- 3 

200 
; 10 

400 



V - 

V - 
V 
V 



K°p+ 

■ir°p+ 
> 7r°e+ 



- 10 

- 1 
0.05 

- 1 



~ 0.01 
~ 1 

~ 0.05 
~ 1 



r(p -> 7TZ/ e ) 



327rm2 w 



v / / 




'x 



Cl2(L e ) 2 l + C 2 \{L^) 12 + (^22(^)21(^)12 



(71) 



where m t = 2.4 • 180 GeV (an artifact of the Al definition). 

Crucial for the suppression of the p — > Kv decay mode is the form of matrix A. We 
had study the renormalization of the matrix A from GUT to the SUSY breaking scale 
and it appears that its form is not changed, so the results presented in this section are 
valid also for the case when the renormalization effects are taken into account. 

Proton decays p — > Kl + and p — > 7re + in our model are described by the box diagrams 
shown on Fig. 2f and 2h and the vertex diagram shown on Fig. 2g amplitude of which 
is proportional to AI/I for the same reason as that described by the diagram of Fig. 2a 
(see the beginning of this section). Vertex diagram shown on Fig. 2e produces b quark, 
so it is irrelevant for proton decay. For proton decay widths we obtain: 



T{p - Kp) 



[m 



l-m\f (I -2a 



327T7B 3 , 



My 



2C 2 \(L e )i 2 + Ci 2 (Lj) 2 l +Ci 2 (L^) 2 ij 

T(p -> Ke) 




m t \V ub \ 



Im 



— (2C 22 + V US C 12 + 
m t 



V ts (c 12 + C 22 (L+) l2 )\ 2 , (72) 



[mt 



2 \2 

k) 



327rm 3 



1 - 2a 



My 



m 2 t \V uh \ 



^ (C 21+ 
m t 



C 22 (L e ) 2 i) — ^— — 



V ts (C 12 (L e ) 21 + C 21 (L+) 12 + C 22 (L+) 12 (L e ) 21 )| 2 , (73) 



T{p -> 7T/i) 




QAirm 



m 2 t \V, 



ah I 



'x 



(2C 12 + C 22 U cd ) - (^J U td (C 12 + C 22 (L+) 12 



/A/ 



(74) 
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Table 3: Numerical values of parameters which were used in the estimates of the proton 
partial lifetimes. 



m w 


m q 


m u 




m s 


m c 


m t 


500 GeV 


500 GeV 


3.6 MeV 


6 MeV 


151 MeV 


1.3 GeV 


180 GeV 


v ub 


V ts 


v td 


sin 2(3 


m) 




a 3 (M z ) 


0.003 


0.05 


0.01 


0.51 


1.9 


10 16 GeV 


0.120 



ne) 




64-Km 



V ub \ 2 x 



m t 



d )12 



C22(L d ~)l2(L e )2i) - 



AI 



V td (Ci2(L e ) 2 l + C2l(L^)l2 + C*22(^)l2(-^e)21 



(75) 



The desirable at GUT scale relations -A = 3 or —3 in our model occur for C02 



61 = 

or C22 = 2A S . In the first case the strong suppression of the p — > Kfi mode will occur. In 
numerical estimates we consider both these cases. 

Proton partial lifetimes for C22 = 2A S and C22 = and for the values of the parameters 
from Table |3] are presented in Table |2| . As we see proton partial lifetimes with emission 



of neutrino for values (rriq, m^) = (500 GeV, 500 GeV) in both cases C22 = 2A S and 
C22 = are compatible with the experimental data. For the case C22 = 2A S the decays 
p — > Kfi and p — >• 7r/x are too fast and we have to change the masses of SUSY particles. 
For example for (m,, m^) = (1 TeV, 100 GeV) we get r(p — > Kfi) = 10 32 years and 
r(p — > 7r/i) = 10 33 years. For the case C22 = the p — > 7r/i mode dominates and for 



[rriq, m w ) 



(600 GeV, 100 GeV) r(p -> nfi) = 10 32 years 



At the end of this section let us present the results of the calculation of the matrix 
elements which contribute into proton decay in our model: 



(uK\(us){ud)\p) 



(vK\{ud){vs)\p) 



[3G 



2V2, 



a 



(M A + M s )/2 3 



(3G 
Ma 



V2{l-2/3a)(vP L p)K 



(vir\(ud)(dv)\p) = ^-V2(uP L p)7r 



(lK\(us)(ul)\p) = ^-V2(l - 2a)(lP L p)K 
(ln\(ud)(ul)\p) = ^(IPlp)* • 



(76) 

(77) 
(78) 
(79) 
(80) 
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6 Discussions 



One of the most appealing next step after minimal standard SU(3) x SU{2) x U(l) 
model is Grand Unification Theory. Both theoretical argument (hierarchy problem) and 
experimental measurements (values of electroweak mixing angle and a s ) prefer, select 
or point out on SUSY GUT. Simplest variant is SU(5) SUSY GUT. However, minimal 
version of the model has two drawbacks: too short proton life time and famous ratio: 
md/ m s — m e / m Li- Proton decay proceed through operators with d — 5. Dominant decay 
mode is p — > Kv with life time of the order of 10 28 years (compare with experimental 
bound r(p — > Kv) > 10 32 years). From the ratio of electron and muon masses we get 
m s /md = 200 which contradicts phenomenological value m s /md = 20 4- 30. 

Both these dissapointing results follow from one source: Yukawa interactions of quark- 
lepton (super) multiplets with Higgs fields in the minimal SU(5). Beyond minimal model 
Yukawa interactions are less constrained. In our approach pattern of quark and lepton 
masses and CKM matrix is explained by the higher dimension operators through which 
first two fermion generations get their masses. Now predictions for proton life time differ 
drastically from that of minimal SUSY SU(5) GUT. Since only third generation fermions 
interacts with 5 and 5 Higgs fields in the same way as in minimal model, operators with 
d = 5 involve these heavy particles which can not participate in proton decay. Bare 
third generation particles get admixtures from first two generations which are small. This 
smallness guarantee smallness of the deviation of CKM matrix from unity. In this way 
proton decay is also suppressed. For scalar quark masses rriq = 500 GeV we obtain: 
r(p — > Kv) ~ 10 32 years which is 4 orders of magnitude better than in minimal model. It 
is interesting to note that p — > Kv decay proceed due to large mass of top quark which 
manifest itself in noticeable mass difference between t and u, c. To suppress proton decay 
further two possibilities exists. First, straightforward one uses heavier squarks. In this 
way p — * K\i mode dominates over additionally suppressed (~ m 2 /m~) p — > Kv mode 
and for rriq = 1.2 TeV and = 100 GeV we get r(p — > Kix) = 10 32 years. Second 
possibility is intimately connected with desirable ratio m^/m s ph ±3 at GUT scale. There 
are two possibilities to get this ratio in our model: to form the 45-plet from the product 
of Higgs fields 24 x 5 (famous Georgi-Jarlskog construction) or to compose 45-plet and 
5-plet in a special way. In the first case contribution to d — 5 operator is of the order 
of \ s while in the second case it is suppressed. It equals zero for m^/m s = 3 and is less 
then 0.1A S for experimentally acceptable choice m M /m s = 2.6 — 3.4. In this way even for 
m 3 = 600 GeV and = 100 GeV we get r p = 10 32 years. 

Search for proton decay at Superkamiokande detector should define future fate of the 
suggested scenario. 
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Figure 1: This picture represents Feynman diagrams which contribute into the leading in 
minimal SUSY SU(5) decay mode p — > K + u fl . Sum of the contribution of Fig. lb and 
Fig. lc equals zero. 
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